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ABSTRACT 



Aims. In this work we investigate the dynamic behavior of inter-network regions of the solar chromosphere. 

Methods. We observed the chromosphere of the quiet Sun using a narrow-band Lyot filter centered at the Ca ii K2v emission peak with 
a bandpass of 0.3 A. We achieved a spatial resolution of on average Cf.'l at a cadence of 10 s. 

Results. In the inter-network we find a mesh-like pattern that features bright grains at the vertices. The pattern has a typical spatial 
scale of r.'95 and a mean evolution time scale of 53 s with a standard deviation of 10 s. A comparison of our results with a recent three- 
dimensional radiation hydrodynamical model implies that the observed pattern is of chromospheric origin. The measured time scales 
are not compatible with those of reversed granulation in the photosphere although the appearance is similar. A direct comparison 
between network and inter-network structure shows that their typical time scales differ by at least a factor of two. 
Conclusions. The existence of a rapidly evolving small-scale pattern in the inter-network regions supports the picture of the lower 
chromosphere as a highly dynamical and intermittent phenomenon. 
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1. Introduction 

Deducing the thermal structure of the solar chromosphere is 
complicated by the limited number of diagnostics and the dif- 
ficulties to interpret them. An example are the core regions of 
the calcium resonance lines (CaiiH and K) that are commonly 
assumed to originate from the chromosphere (cf. lVernazza et al.) 
Il981h . Unfortunately the simplifying assumption of local ther- 
modynamic equilibrium (LTE) is not valid for the emission re- 
versal peaks and the cores of these lines. Rather, non-local ef- 
fects must be taken into account for a realistic description of 
chromospheric plasma and radiation field. This is not only a 
problem for the interpretation of observations of these spec- 
tral features but also for the numerical modeling of the chro- 
mosphere. Th e recent three-dimensiona l radiation hydrodynam- 
ics model by IWedemever et al.l ( l2004 hereafter W04) shows 
a mesh-like pattern at chromospheric heights in non-magnetic 
inter-network regions, consisting of cool regions and hot fil- 
aments that are a product of propagation and interaction of 
oblique shock waves and adiabatic expansion of post-shock re- 
gions. The authors argue that the amplitude of the temperature 
variations in the model chromosphere is somewhat uncertain due 
to the too simple radiative transfer (grey, LTE). But the pattern 
is nevertheless representative as it is a product of acoustic waves 
that are excited in the lower well-modeled layers. The spatial 
scales of the pattern are comparable to the underlying granu- 
lation but the dynamical time scales are much shorter Mapping 
the chromospheric small-scale structure of inter-network regions 
thus requires both high spatial and temporal resolution. This 
might be the reason why most earlier observations failed to de- 
tect this phenomenon (Table 0. 

Here we present new observations taken with a Lyot filter 
of very narrow transmitted wavelength range and a spatial res- 



olution high enough to resolve the chromospheric small-scale 
pattern of inter-network regions. 

In Sect. |2] we describe the observations and data reduction, 
followed by results and discussion in Sects. |3and|3] respectively. 



2. Observation and data reduction 

A high-resolution image sequence was obtained with the 
German Vacuum Tower Telescope (VTT) at the Observatorio 
del Teide on Tenerife, Spain, during UT 8:25-10:21 on April 18, 
2005. The telescope was pointed at a quiet area near disk-center. 
The field of view is 49"x 49" showing chromospheric network 
and inter-network. The data set consists of images observed in 
the G-Band (4305 A) using an interference filter with a FWHM 
of lOA, and, using Lyot filters, in CaiiK (3933A) and Ha 
(6563 A) with FWHMs of 0.3 A and 0.25 A, taken simultane- 
ously. In this paper we only consider the results obtained from 
the calcium filtergrams. A sample image together with closeups 
of a network and an inter-network region is shown in Fig. ^ 
The filter curve of the Ca filter is displayed in Fig. |2] along 



Table 1. Recent observations of the Can resonance lines. The 
properties of the data show that the presented data set is unique 
regarding the combination of spectral and spatial resolution (AA 
is FWHM). 
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Fig. 1. (A) Calcium image with superimposed network mask and closeup regions (numbered squares), (B) the closeups of the 
network region and (C) the inter-network region, (D) the MDI magnetogram, (E) a 2-D power map of the temporal frequency band 
at 5.5 mHz with a width of 0.9 mHz and (F) the map of 1/e time scales. All axes are in arcseconds. 



with filter curves of other observations for comparison. The im- 
age sequence was taken at a cadence of 10 s and a pixel scale 
of O'.'1485/px. There was a need for telescope adjustments at 
UT 9:00, separating the sequence into two parts with a gap of 
2 min. The transmission of the Lyot filter is rather poor, thus the 
light level of the CaiiK line core - being by itself less than 5% 
of the continuum intensity - was reduced even further, leading 
to an exposure time of 2 s. 

All images were corrected for hot pixels and calibrated us- 
ing the standard procedure of flat-fielding. Additionally, the se- 
quence was corrected for the change in mean intensity in each 
image due to the increasing solar elevation with a fitted sine 
curve and scaled to the mean intensity of the time series. Finally, 
the noise in each image was reduced using a Wiener Filter From 
the cutoff of the azimuthally integrated spatial power spectrum 
we estimate a spatial resolution of on average Ql!l or better which 
was achieved by the fair to good seeing (ro ~13 cm) in combina- 
tion with the Kiepenheuer Adaptive Optics System (KAOS). 



3. Analysis and results 

In our analysis we focus on the second half of the image se- 
ries with a length of 75 min. Network regions embody chains 
of persistent bright points (Fig.^J). Inter-network regions show 
a highly dynamic behavior with small-scale, cell-like structures 
(Fig. ^T) and small brightenings at the vertic es that resemble 
bright points or grains (see, e.g., iLites et al]fl999i) . We com- 
puted the spatial power spectrum of the inter-network pattern. 
In order to enhance the small-scale nature of the pattern of in- 
terest, an unsharp-masked image was subtracted from the data. 
The two-dimensional power spectrum of the resulting image was 
then azimuthally integrated. The one-dimensional power spec- 



trum is shown in Fig. The barycenter of power is at spatial 
frequency //, - 0.512 arcs"' corresponding to r.'95. The bright 
grains have typically a size of 1-2 "(see Fig.^^)- 

To demonstrate the dynamical behavior in further detail, we 
calculated the two-dimensional autocorrelation function for sub- 
fields of r.'485 X r.'485 which overlap by half of their field size 
in both directions. The evolution time scale is defined as the 
time at which the autocorrelation has dropped to a value of 1/e. 
Finding the evolution time scale for each subfield separately re- 
sults in an array of 66 x 66 values over the field of view. The map 
in Fig.^ shows the results for the individual subfields. Figure|3 
shows a plot of time scale versus mean intensity for all subfields. 
In order to persue a more detailed analysis of the different behav- 
ior of network and inter-network, criteria need to be established 
which distinguish these two regions. 

As a first step, the lower mean brightness of the inter- 
network pattern was used to create a mask of network and inter- 
network regions b y thresholding the smoothed temporal mean 
of the intensity (cf.'Ru tten et alJ2004h . We found such a thresh- 
old at 1=1.017 by fitting the sum of two Gaussians to the inten- 
sity histogram of the mean image. The threshold is indicated as 
vertical line in Fig. |3 which divides the plot into two regions, 
the left containing dark subfields, the right containing the bright 
subfields. The application of the threshold in combination with 
a morphological operation to remove isolated pixels resulted in 
mask Ml which is indicated in Fig.^^, D, E and F as a solid 
contour line. Interestingly, the time scale distribution of sub- 
fields in Fig. 13 shows a gap between 100 s and 200 s that sep- 
arates long and short time scales (shaded area). We therefore 
chose a threshold of 150 s to discriminate between the fast and 
the slow evolving regimes (horizontal line in Fig. Together 
with the intensity threshold the new time scale criterion divides 
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Fig. 2. The calcium line around 3933.5 A llNeckellll999l) (thin 
sohd) plotted together with theoretical Lyot filter transmission 
curves with 0.3 A (thick solid), 0.6 A (dashed) and a Gaussian 
with 1.350 A (dot-dashed) FWHM (see Table[D. 



the plot into quadrants. Applying the time scale threshold to the 
image sequence results in a new mask M2, which is indicated 
as dot-dashed contour in Fig.^. Mask M2 is almost completely 
contained within mask Ml. The pixels that are attributed to the 
network in mask Ml but not in M2 belong to subfields which 
contribute to the lower right quadrant of Fig. ^ This quadrant 
contains subfields that are bright on average and exhibit short 
time scales. Apparently, there are many structures which belong 
to the inter-network even though they are brighter than average. 
That includes subfields with bright CaiiK grains, implying a 
connection between grains and inter-network pattern. The ma- 
jority of fast evolving subfields, however, is darker than average 
and clearly belongs to inter-network regions (lower left quad- 
rant). The barycenter of these subfields (with time scales be- 
low 150 s) Ues at = 53 s and a normaHzed intensity of 0.953. 
The standard deviation in time scale is 10 s, the minimum value 
is ^pe.min - 39 s. In Contrast, network structures show long time 
scales, some of which were longer than the duration of our data 
sequence, and are on average brighter than inter-network struc- 
tures. Therefore, they are found mostly in the upper right quad- 
rant of Fig.0] As expected, these persistent patches of enhanced 
intensity match magnetic structures shown in the MDI magne- 
togram (Fig. ^3) identifying them as part of the network. There 
are only few dark and slowly changing structures which appear 
in the upper left quadrant. 

We computed temporal power spectra for each pixel and av- 
eraged over network and inter-network pixels according to mask 
M2. The result in Fig.|5lshows a clearly enhanced signal in the 
frequency band around 5.5 mHz within inter-network regions. 
This can also be seen in the power map of Fig. ^ where the 
band around 5.5 mHz with a band pass of 0.9 mHz is displayed. 
The network mask is again outlined by contour lines. The pixels 
of decreased power are almost completely contained within M2. 



4. Discussion and Conclusion 

We provide a new method to discriminate network and inter- 
network structures in filtergrams which show chromospheric 
structure when the data are not accompanied by information 
about the underlying magnetic field. Figure |5 demonstrates that 
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Fig. 3. One-dimensional power spectrum of the inter-network 
pattern. 



the pattern evolution time scale provides an excellent criterion to 
distinguish between network and inter-network regions: all inter- 
network points are grouped in the time interval below 100 s, 
whereas the network has time scales above 200 s. Temporal 
thresholding delivers an accurate mask, M2, for network pixels. 

Structures in the inter-network similar to those found in 
our observations have already b een detected in various other 
works. Leenaarts & We demever- Bohm (2005, hereafter LW05) 
investigate a similar pattern at higher spatial resolution - but 
lower spectral resolution - and identified it with reversed gran- 
ulation. This result is also supported by the much longer 
time scale of 2min they determine for the pattern. In Fig. 1 
of Kneer & von Uexkiill ( 1993) the pattern can be surmised. 
Tritschler & Schmidt (2002) present observational data with 
time scales of inter-network structures close to the ones obtained 
in this work. Their spatial resolution did not suffice to resolve the 
small-scale pattern and they used a broader bandwidth of 0.6 A. 
This work is the first one to observe the small-scale pattern with 
a mean evolution time scale much shorter than 100 s. 

In order to verify the chromospheric origin of the detected 
inter-network pattern, we compare its spatial and evolution time 
scale to the recent 3D model by W04. The temperature ampli- 
tudes of the model chromosphere might be somewhat uncer- 
tain due to the simplified treatment (grey, LTE), whereas the 
excitation of the waves occurs in realistically modelled lower 
layers, making the dynamics and related time scales more re- 
liable. A direct comparison with synthetic intensity maps for 
Ca n K2v is not possible yet as the model chromosphere is dom- 
inated by strong shocks that cause numerical problems. We can, 
however, perform a qualitative comparison with horizontal slices 
of gas temperature from the model chromosphere. The cuts ex- 
hibit a pattern of granular size at a geometrical height around 
1000 km above r = 1 which is rapidly evolving with time scales 
of around = 20 - 30 s. The pattern is quite similar to the one 
we observe, but the evolution time scale is still a factor of two 
shorter than the observational value. This can be explained as 
follows: First of all, the correlation between gas temperature at 
z = 1000 km and intensity in the CaiiK2v spectral feature might 
be poor as the assumption of LTE is not valid in this strongly 
scattering line. In LTE local changes in temperature are instanta- 
neously translated into changes of the source function and thus 
intensity changes although the intensity response is non-linear 
in amplitude. In contrast to the temperature slices at a well- 
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Fig. 4. Time scale versus intensity. The gap between 100 and 
200 s indicated with a shaded horizontal band shows the different 
behavior of network and inter-network structures. 



defined geometrical height, the emergent intensity (at a given 
wavelength) is a non-local quantity that has contributions from 
an extended height range. For the line core, the main contribu- 
tion of intensity originates from chromospheric height s around 
1000- 1500km above t=1 (cf. IVemazza etalJll981h . but al- 
ready at 0.15 A off-core considerable contributions to the filter- 
gram intensity from lower geometrical heights have to be ex- 
pected (Rammacher, priv. comm.). The emergent intensity is 
thus integrated along the line of sight over many layers which 
individually would sho w different pattern evolution times (de- 
creasing with height, see Wedemever 2003). The integrated pat- 
tern is consequently smoothed so that the relative intensity vari- 
ations become smaller and result in longer evolution time scales. 
This effect is even enhanced due to the fact that the observed in- 
tensity refers to corrugated surfaces of optical depth instead of 
plane horizontal cuts as used for the model. The intensity pattern 
will be smoothed even more compared to the fine structure in 
the temperature cuts. W04 calculate a time scale of ^pe,z > 120 s 
for changes in gas temperature in a plane horizontal cut at the 
bottom of the photosphere, whereas the grey emergent intensity 
(under assumption of LTE, which reflects conditions at the bot- 
tom of the photosphere), yields an increase of the typical time 
scale to ^pe,T = 200 s. This is due to the fact that variations 
tend to be smaller on surfaces of optical depth because these are 
shaped by inhomogeneities themselves. In addition, seeing and 
instrumental effects lead to a spatial smearing of the intensity 
signal and thus to longer time scales. We found that smoothing 
the temperature cuts to the spatial resolution encountered in our 
data increases the time scales by about 25%. The simulations 
of W04 do not include magnetic field. The existence of inter- 
network fields (with intrinsic field strengths varying from 1 kG 
down to ^lOOG) is undisputed as they hav e been recently mea- 
sured with high polarimetric precisi on (see ISocas-Navarro et aP 
120041: lOtes & Socas-Navarr o'2004\ and references therein). The 
contribution of possibly unresolved magnetic elements in inter- 
network regions tends to lengthen the effective evolution time 
scale. A preliminary analysis of temperature cuts from the 3D 
magnetohydrodynamics model by Schaffenberger et al. (2005) 
results in a slightly higher value of fpe.MHD = 27 s at the same 
height (z - 1000 km) compared to 24 s in the purely hydrody- 
namical model by W04. Consideration of all these effects leads 
to a reasonable agreement between observation and numerical 
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Fig. 5. Temporal power for network (dashed) and inter-network 
(solid). 



model. An interpretation of the pattern as reversed granulation 
is inappropriate as the corresponding time scale is much longer 
than the 53 s found here. LW05 state 120 s and 90 s for observed 
and synthetic CaiiH line wing intensities, respectively, which 
clearly originate from the mid-photosphere (LW05, Wedemev^ 
2003). These values are in line with the time scale of order of 
90 s for the reversed granulation pattern observed in the Fe I line 
core at 709.0 nm by Janssen & Cauzzi ( 2006). We conclude that 
the small-scale inter-network pattern reported here evolves too 
fast to be due to reversed granulation and is thus of chromo- 
spheric origin. The pattern is very likely the intensity signature 
of the mesh-like, fast evolving pattern found in recent 3D hydro- 
dynamical models of the solar chromosphere. It is presumably 
a product of propagation and interaction of shock waves (W04) 
that are excited below th e observed layer The sa me is true for 
the Can grains (see also ICarlsson & Steinlll99'7h that we con- 
sider an integral part of the pattern. Collision of two neighbour- 
ing wave fronts results in compression and heating of the gas 
at the contact region which shows up as a band of enhanced in- 
tensity, making up the "strings" of the mesh-like pattern. At a 
mesh vertex, however, more than two wave fronts meet, result- 
ing in stronger compression and thus in an intensity higher than 
the remaining pattern. Due to the limitations of earlier observa- 
tions (see Tabled only the brightest components have been de- 
tected as grains whereas the pattern itself remained undetected. 
A more detailed study including the simultaneous G-band data 
is in preparation. 
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